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The Dalitz plot slope parameters g, h and k for the K decay have been measured 

using in-flight decays detected with the ISTRA+ setup operating in the 25 GeV negative sec- 
ondary beam of the U-70 PS. About 252 K events with four- momenta measured for the tt~ and 
four involved photons were used for the analysis. The values obtained 

g = 0.627 ± 0.004(stat) ± O.OlO(sysi), 
h = 0.046 ± 0.004(stat) ± 0M2(syst), 
k = 0.001 ± O.OOl(siai) ± 0.002(syst) 



are consistent with the world averages dominated by K + data, but have significantly smaller 
errors. 



1 Introduction 



The determination of the Dalitz plot slope parameters for the K — > (Sti)^ decays is 
of interest for the Chiral Perturbation Theory (ChPT) and as a check on the direct CP 
violation. The latter would manifest itself by the difference between the K + and K~ 
Dalitz plot distributions (see for example pQ). 

The square of the matrix element of the r' (K — *> ir ^^) decay can be written as 

{Atf* -> 3tt)| 2 oc 1 + gY + hY 2 + k X 2 + ... , (1) 

where X = (si—S2)/m 2 [ andF = (53— so)/m 2 T are the Dalitz variables, and the parameters 
g k are the "Dalitz plot slopes". Here S{ = (px — Pi) 2 , so = ~(si + s 2 + s 3 ), px and Pi 
are the K ± and 7Tj four- momenta (7r 3 is the odd pion). 

The direct CP violation could be detected by the observation of the following charge 
asymmetry: 

ifig)r> = 9 4^-- (2) 

9 + + 9 

The theoretical predictions for the asymmetry (Sg) T i in the framework of the Standard 
Model (SM) were originally spread in the wide range of ~ 2 ■ 10~ 6 -v- 10 -3 |2j- Over 
last years they have converged to the value of ~ 10~ 5 [HE]. In a wide class of possible 
super symmetric extensions of the SM larger values are possible. For example, in the 
Weinberg model jlj the value of ~ 2 ■ 10~ 4 is predicted 0. 

This topic has been attracting significant interest in the recent years. The comparison 
of the results of the latest K + experiment [0] with the only one existing K~ measurement 
[7j leads to the value of (Sg) T i = 0.117 ± 0.022, i.e. 5.3 sigma effect. 

Recently the CP conserving amplitudes for the K — ► 37r decays were recalculated in 
ChPT at the next-to-leading order [S]. In this view a new high statistics measurement 
for the K — > 3ir decays is desirable. 

These observations encourage us to perform a new measurement of the Dalitz plot 
slope parameters g, h and k for the K~ decay, based on the statistics of about 

250 K events. 



2 Experimental setup 

The experiment has been performed at the IHEP proton synchrotron U-70 with the ex- 
perimental apparatus "ISTRA+" which is a modification of the "ISTRA-M" setup 
and described in some details in our recent papers [TU], where studies of the K~ 3 and 
K~ 3 decays have been presented. The setup is located in the negative unseparated sec- 
ondary beam with the following parameters during the measurements: the momentum is 
~ 25 GeV/c with Ap/p ~ 2 %, the admixture of kaons is ~ 3 %, and the total intensity 
is ~ 3 ■ 10 6 per spill. 

The side view of the "ISTRA+" detector is shown in Fig.QJ The momentum and the 
direction of the beam particle, deflected by the magnet Ml, is measured with four pro- 
portional chambers BPC1-BPC4. The kaon identification is done by three threshold gas 
Cherenkov counters C0-C2. The momenta of the secondary charged particles, deflected 
in the vertical plane by the magnet M2, are measured with three proportional chambers 
PC1-PC3, with three drift chambers DC1-DC3, and with four planes of drift tubes DT. 
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Figure 1: The side view of the "ISTRA+" detector. 



The secondary photons are registered by the lead glass electromagnetic calorimeters SP1 
and SP2. To veto low energy photons the decay volume is surrounded by the guard 
system of eight lead glass rings and by the SP2. The wide aperture threshold helium 
Cherenkov counters C3 and C4 are used to trigger the electrons and are not relevant for 
the present paper. In Fig.Q HC is a scintillator-iron sampling hadron calorimeter, MH is 
a scintillator hodoscope used to solve the reconstruction ambiguity for multitrack events 
and improve the time resolution of the tracking system, MuH is a muon hodoscope. 

The trigger is provided by the scintillation counters S1-S5, the Cherenkov counters 
C0-C2, and the sum of the amplitudes from the last dinodes of the calorimeter SP1 (see 
Ref. jTU] for details). The latter serves to suppress the K~ — > i^v^ decay. 

3 Event selection 

About 363 M and 332 M events were collected during two physics runs in Spring 2001 and 
Autumn 2001. These experimental data are supported by about 260 M events generated 
with the Monte Carlo program GEANT3 [11J. The Monte Carlo simulation includes a 
realistic description of the experimental setup: the decay volume entrance windows, the 
track chamber windows, gas, sense wires and cathode structures, the Cherenkov counter 
mirrors and gas, the showers development in the electromagnetic calorimeters, etc. The 
details of the reconstruction procedure have been presented in Ref. ^U]; here only key 
points relevant to the K~ — > 7r~7r°7r° event selection are described. 

The data processing starts with the beam particle reconstruction in the proportional 
chambers BPC1-BPC4, and then with the secondary tracks reconstruction in the tracking 
system PC1-PC3, DC1-DC3 and DT. Finally, the electromagnetic showers are looked for 
in the calorimeters SP1 and SP2. The method of the photons reconstruction is based on 
the Monte Carlo generated patterns of showers. To suppress leptonic decays of kaons the 
particle identification is used JUj- The muons are identified using the information from 
the calorimeters SP1 and HC. The electrons are identified using the ratio of the energy 
of the shower, detected in the SP1 and associated with the track of the electron, and the 
momentum of the electron. 

At the first step of the event selection only the measurements of the beam and sec- 
ondary charged particles are used. Those events are selected which satisfy the following 
requirements: 

- only one beam track and only one negative secondary track are detected; 

- the probability of the vertex fit, CL(x 2 ), is more than 0.01; 



- the decay vertex is in the region before the calorimeter SP2 (400 cm < z < 1650 
cm) , and its transverse deviation from the setup axis is less than 10 cm; 

- the angle between the beam track and the secondary track is more than 2.5 mrad; 

- the transverse momentum of the secondary track with respect to the beam direction 
is less than 150 MeV/c; 

- the secondary track is not identified as an electron or as a muon. 

At the second step of the event selection the measurements of the showers in the 
calorimeters SP1 and SP2 are used. The following requirements are used to choose the 
photons (showers): 

- the distance between the shower in the calorimeter SP1 and the intersection of the 
secondary track with the transverse plane of the SP1 is more than 1.5 cm along the 
magnetic field of the M2 and more than 9 cm in the transverse plane; 

- the photon energy is more than 0.7 GeV, but is more than 1.4 GeV when the photon 
is detected in three or less cells of the calorimeter SP1; 

- for events where the secondary track is not associated with any shower in the 
calorimeter SP1 and with any hit in the hodoscope MH the photon energy is more 
than 1.4 GeV; 

- the energy of the photon found in the multishower cluster is more than 2 GeV. 

For each selected 77 pair the deviation of its effective mass from the 7r° mass, AM(77) = 
M(77) — m(7r°), is calculated. If this deviation is in the range of |AM(77)| < 50 MeV, 
the 77 pair is considered as a candidate of the 7r° decay. Then, if this 77 pair is taken as a 
7T° decay, the four-momenta of its photons are multiplied by a factor A = m(7r°) / M (jj) . 
The 7T° detection is illustrated in Fig.Et, where the spectrum of the effective masses of 
the "best" 77 pairs, i.e. the pairs with the smallest value of |AM(77)|, is shown for the 
selected events with at least four detected photons. In further selections all combinatorial 
77° candidates are used, not only the "best" . 

At the third step of the event selection the sample of the K~ — > 71-^7777 events with 
two 7T° — > 77 decays is collected. For this sample the further selection is done by the 
requirements that the measured value of the kaon mass is in the range of |M(7r _ 7r 7r°) — 
m(K~)\ < 80 MeV (see Fig.|2b) and then the event passes the kinematical 6C-fit for the 
K~ — > 7r~7r°7r° hypothesis. The efficiency of the last cut is about 89%. Both requirements 
are considered for all possible combinations of photons and the best 6C-fit hypothesis is 
chosen. 

Using the mentioned above selection criteria for the K~ — > 7r~7r 7r° decay we have 
collected 252 K completely reconstructed events. The corresponding numbers of accepted 
Monte Carlo events are about ten times larger than in the experiment. The surviving 
background, which is mainly due to kaon decays to other allowed modes, is estimated 
from the Monte Carlo simulation to be less than 0.04%. 

The detailed event reduction statistics is given in Table ^ The difference in the 
fraction of the selected events for two runs is due to the higher threshold in the SP1 at 
the trigger level for the second run, i.e. the initial statistics of the first run contains more 
K~ — > i^v^ decays than the initial statistics of the second run. 
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Figure 2: a) The deviation AM(77) = M^yy) — m(ir°) of the effective mass of the 77 pair 
with the smallest value of |AM(77)| in the selected events with at least four detected photons, 
b) The deviation AM(7r~7r°7r°) = M(7r~7r°7r°) - m(K~) of the effective mass of the 7r"7r 7r 
system with the best photon combination in the selected events. 



Table 1: The event reduction statistics for two physics runs. 



Run 


Spring 2001 


Autumn 2001 


Total number of events 


363 M 


332 M 


Beam track reconstructed 


269 M 


248 M 


Secondary track(s) reconstructed 


121 M 


124 M 


Number of events written on DST 


93 M 


108 M 


K~ and tt~ selected (all cuts of the 1st step) 


4882 K 


7432 K 


at least four 7 selected (all cuts of the 2nd step) 


83 K 


244 K 


two 7T° selected 


69 K 


218 K 


K~ — > 7r~7r°7r° selected (all cuts of the 3rd step) 


59 K 


193 K 
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4 Analysis 



To determine the parameters g, h and k in Eq. (jlj the uncorrected distribution p(X, Y) of 
the event density on the Dalitz plot was analysed. This distribution is shown in Fig. El for 
the selected K~ — > 7r~7r°7r° events. At first the background contamination was subtracted 
from the Dalitz plot. The contamination was estimated from the Monte Carlo simulation 
of the particle interaction with the material of the detector and of the kaon decay including 
all decay modes with the branching ratios more than 1%. The corresponding branching 
ratios and matrix elements in this simulation were taken from the PDG [T2]. In Fig.0] 
the fraction of the background contamination is shown as a function of the Dalitz plot 
variables IX I and Y. 
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Figure 3: The Dalitz plot, Y = (s 3 
K~ — > 7r~7r°7r° events (uncorrected). 



so)/m^ versus X = (si — S2)/m^, for the selected 




Figure 4: The fraction of the background contamination as a function of the Dalitz plot variable 
estimated from Monte Carlo simulation. 

The background subtracted distribution p'(X,Y) was fitted by the method of least 
squares with the function: 

p'(X,Y) oc F 1 (X,Y)+gF 2 (X,Y) + hF 3 (X,Y) + kF 4 (X,Y), (3) 
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where Fk(X,Y) are the distributions for the w^— weighted Monte Carlo K~ — > 7r~7r°7r° 
events generated with the constant matrix element and reconstructed with the same 
program as for the real events. The weight factors w\ = 1, w-i = Y true , u> 3 = Y^ rue 
and w 4 = Xf rue are given by the "true" values of X and Y, but the bins of F k (X, Y) are 
given by the "measured" ones. This method allows to avoid the systematic errors [T5] 
due to the "migration" of the events on the Dalitz plot because of the finite experimental 
resolution. Fig.|3] illustrates the Monte Carlo estimated experimental resolution for the 
variables \X\ and Y. where the "measured" values are shown versus the "true" ones. 
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Figure 5: The "measured" values of the Dalitz plot variables \X\ and Y versus the "true" ones 
estimated from Monte Carlo simulation. 



5 Results 

The result of the least squares fit is illustrated in Fig-El where the matrix element 



\A(K-^n-n\°)\> = C • (4) 

for the K~ — > 7r~7r 7r° decay is shown as a function of Y in the different intervals of |X|. 
The normalization constant C in Eq. (jlj) provides the value of \A\ 2 = 1 at the point (X=0, 
Y—0). The integrated dependences of the matrix element on the variables Y and |X| are 
shown in Fig.0 

The values of the Dalitz plot slope parameters are found to be 

g = 0.627 ±0.004 ±0.010, 
h = 0.046 ±0.004 ±0.012, 
k = 0.001 ±0.001 ±0.002 

with x 2 / n df = 502/558. Here the first errors are statistical and the second ones are 
systematic. 

In the determination of the systematic errors of the slope parameters measurement 
the following contributions were taken into account. 
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Figure 6: The matrix element dependence on the variable Y in the different intervals of the 
variable |X| for the K~ — ► 7r~7r°7r° decay. The curves are the result of the fit to the function 

©■ 
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Figure 7: The integrated dependences of the matrix element on the variables a) Y and b) |X| 
for the K~ — > 7r _ 7r 7r° decay. The curves are the result of the fit to the function (j5J). 
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- Two samples of the events collected in the runs (Spring and Autumn 2001) with 
some differences in characteristics of the setup were analysed separately (the corre- 
sponding contributions to the systematic errors are Ag = 0.001, Ah = 0.001 and 
Ak = 0.0002). 

- To avoid some uncertainties at the edge of the Dalitz plot the extreme bins of this 
plot were cut (Ag = 0.001, Ah = 0.002, Ak = 0.0004). 

- The energy threshold of the selected photons was increased from the value of 0.7 
GeV to 2 GeV (Ag = 0.006, Ah = 0.007, Ak = 0.001). 

- The mass and energy ranges used in the event selection criteria were varied from the 
value of 30 MeV to 80 MeV with and without introducing a factor A = m(TT°)/M('yy) 
(Ag = 0.003, Ah = 0.004, A A; = 0.001). 

- The background contamination estimated from the Monte Carlo simulation was not 
subtracted from the Dalitz plot before the least squares fit (Ag = 0.001, Ah = 
0.001, Ak = 0.0002). 

- The particle identification of the secondary track was not used (Ag = 0.001, Ah = 
0.001, Ak = 0.0003). 

- The electromagnetic showers in the calorimeter SP2 were not used in the photon 
reconstruction (Ag = 0.006, Ah = 0.007, A A; = 0.0014). 

- The upper edge of the decay vertex position was varied along the setup axis between 
the chamber PCI and the calorimeter SP2 (Ag = 0.004, Ah = 0.005, Ak = 0.001). 

- The cut applied to the angle between the beam and secondary tracks was varied 
from 1.5 to 3.5 mrad (Ag = 0.001, Ah = 0.001, Ak = 0.0003). 

6 Summary and conclusion 

The Dalitz plot slope parameters for the K~ — > 7r~7r°7r° decay have been measured using 
the "ISTRA+" spectrometer. It is an update of our preliminary presentation [Tl]. The 
results of our measurement, the world averages [12] and the results of previous experiments 
[6, 7, 15 - 21] on the K ± — > -k ± -k ti decays are presented in Fig.lHl Among the previous 
experiments there are eight measurements of the K + decay, but only one of the K~ 
decay. Our values of the slope parameters g and h are consistent with the world averages 
dominated by K + measurements. The difference between the values of the linear slope 
g obtained for the K~ — > 7r~7r°7r° decay in our experiment and in another one [7] is 1.9 
standard deviations. 

Using the same rules as in the PDG ^2] (without the data (201 E] obtained from the 
linear fit only) and adding our measurement, one can obtain the world average values of 
the linear slope g for the K + and K~ decays separately: g + = 0.684 ± 0.033 and g~ = 
0.617±0.018. They give for the charge asymmetry (j2J) the value of (Sg) T r = 0.051 ±0.028. 

The analytical expressions of the CP conserving K —>■ 3n amplitudes with single 
parameter functions, recalculated in the framework of ChPT at the next-to-leading order, 
were fitted in Ref. [B] to all available K — > 2-k and K — > 3ir data. The result of this global 
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Figure 8: The Dalitz plot slope parameters g and h for the K — > ir 7r°7r° (solid circles), 
K + — ► 7T + 7r 7r (open circles) and K ± — > 7r ± 7r 7r° (solid stars) decays. 

fit gives the following predictions for physical observables of the K ± — > tt^tt ^ decays: 
g = 0.638, h = 0.074 and k = 0.0045. These predictions agree within 1 2 standard 
deviations with our results. 

The work is supported by the RFBR contract no. 03-02-16330. 
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